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We present the simplest model to which one can apply the supersym- 
mctric Hubbard operators recently introduced. [1] For the atomic model, 
H = —EdX 00, where Xoo = |0 >< 0| is a Hubbard operator and Ed is 
the energy of the localized spin level, we show how one can develop exact 
solutions for the entropy and heat capacity as a function of temperature. 
With this gold standard we are able to develop a controlled approximation 
scheme to field theoretically treat the susy approximation at the level of 
mean field + gaussian corrections and test its accuracy against the widely 
used slave boson and slave fermion approximations. We find that in ad- 
dition to slave boson and slave fermion solutions, a new class of solutions 
exists in the physical case Q = 1, N = 2 which can be properly treated 
by neither previously existing approach. The phase diagram generated by 
the mean field saddle-point bears a superficial resemblance to the V-shaped 
phase diagram common to systems close to a quantum critical point and 
may provide a natural starting point for investigations of strongly corre- 
lated models capturing this physics. 

PACS numbers: 71.27.-Ha, 71.10.Hf, 75.20.Hr, 75.40.-s 



1. Introduction 

In heavy fermion systems close to an antiferromagnetic quantum criti- 
cal point (aqcp) one now has the ability to tune from an antiferromagnetic 
state to a paramagnetic state as a function of field [2], pressure [2, 3] or 
doping[2, 3]. In the region above such a point, anomalous behavior is seen 
in measurements of heat capacity ^ ~ ln{^) and resistivity p ^ T. The 
question has arisen in the literature whether these exponents are due to 
a dynamical d=2 local qcp[4] or rather a symptom of a breakdown of the 
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d=3 Fermi liquid[5]. In particular, the observation of hyperscaling in these 
compounds has been used to argue the need for a microscopic theory capa- 
ble of describing the breakdown of the Fermi surface and the formation of 
antifcrromagnctic order as one tunes through this point. 

Susy Hubbard operators[l] with SU(l|l)xU(l) symmetry aUow one to 
tune from a bosonic description of a spin (magnetic) to a fermionic de- 
scription (Fermi surface). For an atomic model, we show using counting 
arguments how it is possible to generate exact solutions for the entropy and 
heat capacity as a function of temperature for the Hilbert space of these 
operators. Comparison of these with field theoretic results (mean field + 
gaussian fluctuations) allows us to estimate the error endemic to this 
approach, contrasting it with slave fermion and slave boson approaches. 

2. An exact solution 

To fix an irreducible representation of the Hubbard operators, we set 

\a >< b\ = Xab = BiBk + F^Fb (1) 

where Fa = (/i, /at, 0) and Ba = (61, ^at, x) defines spin fields to be 
fermionic or bosonic respectively and their slave partners the converse, while 
maintaining the constraints Q = rih + + Uf + n^^ and Y = n^p + Uf — 
(ub + n^) + ^ [0, 0^] where = J2a ^tfa- — is an operator interconverting 
fermions and bosons for the corner state. This generates a series of L-shaped 
Young tableaux, the simplest of which (a single box) corresponds to a single 
physical spin when N=2. To find the exact free energy of the state (Q,Y) 
we simply count the number of available states to the system at a given 
energy level. Defining h = '^^^^^ and w = , the number of states 

available are: 

«.,»,^(-)(---),,^ ,) 

For an atomic model, the Hamiltonian is given by 

H = EdXa-a = EdQ — EdXoo = —E^Xoo (3) 

where E^i is the energy of the d or f-clectron state (we have dropped a 
constant in the free energy) , which leads to a partition function of the form 

Z = J2j2^N{h-i,w- j)e^^i+^^^'i (4) 

j=0 j=0 
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and the free energy F = — rin(Z), entropy S = —BtF and heat capacity 
Cy = T&rS immediately follow. Including constraints, at the mean field 
level the Hamiltonian becomes 

H = - Ed{n<j> + %) + A(n5 + n/ + + - Qq) 

+ C{hf + n^~ (fib + %) + 7^ < [0, e^] > -Yo) (5) 

where we can only evaluate the last term at the level of gaussian fluctuations. 

Nonetheless we include it in the saddle-point, as we additionally treat the 
effects of gaussian fluctuations in the bosonic/fermionic character, 

F = NFf + NFb + F^ + F^ + Fr,- XQo - (Yq + Fgx^ + Fgx, (6) 

which at the mean held level sets 

nf + jfn^ + j^il - Ua) = h (7) 
m + jfUy^ + j^ria = w (8) 

where n„ = hf = ^0(x])^^ , = ^^P^'^^ = ^-f^(Ea\)_i 

riy^ = ^_g(g^^_A^)_|_^ - Analytic solution of these equations leads to the phase 
diagrams shown in Fig.l for the special case N=2. For general I shaped 
Young tableaux one recovers the mean held results of slave fermions (ver- 
tical) and slave fermions (horizontal) although the free energy contains di- 
vergent terms in this limit. 

While the appearance of a mixed phase in the phase diagram of Fig. lb is 
quite suggestive given that the mean held entropies also match along these 
lines, one is ultimately disappointed when one realizes that: i) although wc 
expected that the magnetic phase would require a bosonic description of the 
spins-slave fermions and correspond to more tightly bound spins (Ed<0) 
and slave bosons provide a natural candidate for a heavy Fermi surface 
which might be expected to appear close to (Ed=0) this does not seem 
to be the case here; ii) as the slave fermion mean held over-estimates the 
entropy at N=2, at the mean field level the mixed phase (interpolating to 
the slave boson mean field) ends up having a negative heat capacity. This 
problem appears to be over-come by including the gaussian corrections, 
but in the limits h ^ 1 or w ^ 1 the gaussian fiuctuations are not well- 
defined. If after removal of the non-physical divergences the entropies of the 
mixed phase still match those of slave boson and slave fermion along the 
phase boundary lines (as was suggested by the mean field), then the positive 
difference Ssiavehoson{T oo) -SgiavefermioniT = j^) would imply a smah, 
positive heat capacity for this interval, in accordance with the exact result. 
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Fig. 1. The analytic saddle-point (s-p) solutions at N = 2 allow two possible phase 
diagrams (lines indicate equal free energies): a) The physical s-p picks the smallest 
free energy of the saddle-point solutions, b) A non-physical s-p which picks the 
highest free energy. While minimization of the free energy does choose the correct 
physical result, this second solution illustrates how a non-trivial phase diagram 
may arise from a symmetry of the underlying formalism. In an interacting model 
one might imagine that if the exact solution has ^ ?a In(^) then perhaps in that 
region a similar phase diagram to b) may be the energetically favourable one. 

iii) even were this the case, one would have to accept an ansatz choosing 
the maximal saddle-point free energy to admit Fig. lb. 

In conclusion, wc have shown that for a simple atomic model one re- 
covers the mean field constraints known from slave boson and slave fermion 
treatments. We have shown that a non-trivial mixed solution exists even in 
the physically relevant case Q=l, N=2. Study of the properties of mixed 
solutions in the controllable large N limit may be of interest, and for the 
atomic model can be compared with exact results shown here. A more 
thorough treatment will be given in the near future. [6] 
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